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Abstract. The ppN and ppA three-body systems have been studied within the framework of the fixed 
center approximation of Faddeev equation. The pp interaction in isospin I — 0, spin S = 2 is strongly 
attractive, and so are the Np, Ap interactions. This leads to bound states of both ppN and ppA. We find 
peaks of the modulus squared of the scattering matrix around 2227 MeV for ppN, and 2372 MeV for ppA. 
Yet, the strength of the peak for the ppN amplitude is much smaller than for ppA, weakening the case 
for a ppN bound state, or a dominant ppN component. A discussion is made on how these states can be 
searched for in present programs looking for multimeson final states in different reactions. 

PACS. 14.20. Gk Baryon resonances (S=C=B=0) - 21.45.-v Few-body systems 



1 introduction 

While it has been always accepted that the meson baryon 
cloud on top of three constituent quarks plays some role 
in baryon observables, like form factors and static proper- 
ties, it is becoming gradually more broadly accepted that 
in some cases this meson baryon cloud can even account 
for the largest part of the wave function, as far as one is 
concerned with physical processes taking place at small 
or intermediate energies. The use of effective theories that 
provide interaction Lagrangians in terms of the observable 
degrees of freedom, mesons and baryons [1], has allowed 
to study systems made out of mesons and baryons, and 
in some cases the interaction is so strong that one gets 
bound states or resonances of that system. The systematic 
use of chiral Lagrangians [lll5J|31ll] and unitary techniques 
in coupled channels has lead to the so called chiral uni- 
tary approach which has been rather successful describing 
J^ — l/2~ low lying baryonic resonances ^ or the low 
lying J^ = 3/2^ states [6l[7]. Extensions have been done 
for the interaction of vector mesons with baryons, both 
for the octet of baryons [5] as for the decuplet |2] ■ 

Systems of two mesons and one baryon have also started 
to be studied. Faddeev calculations of these systems have 
found that the low lying J^ — 1/2+ baryonic states, up to 
the Roper which is a very complicated object [TD] , can be 
reproduced as bound states or resonances of two mesons 
and one baryon in coupled channels [TT1I12J . 

States of larger spin have also been investigated, in- 



with positive parity have also been investigated as bound 
states or resonances of two mesons and one baryon. For in- 
stance in Ref. [14j a plausible explanation was found for a 
Z\(l/2+) state around 1750 MeV, which is highly problem- 
atic in quark models [15) . in terms of the wpA system. It 
is interesting to pursue this search, looking for states with 
higher spin which are also difficult to obtain with quark 
models, and which can be easier accommodated by means 
of meson baryon components. In this paper we go along 
this line and study systems made out of two p mesons and 
a. N or A. 

In order to evaluate the amplitudes for the interaction 
of these three body systems we take advantage of the fact 
that the pp interaction in / = and spin S* = 2 is huge 
and leads to the bound state of the /2(1270) resonance, 
which decays into tttt as found in Ref. 16]. The system is 
so bound and packed because of the short range of the in- 
teraction that the presence of an extra baryon should not 
alter very much this strongly bound pp system. This is 
an ideal situation to apply the Fixed Center Approxima- 
tion (EGA) to the Faddeev equations ^17,18,19,20 , which 
certainly simplifies the technical work. Since we are not in- 
terested in a precise determination of the energies, and un- 
certainties of even 50 MeV are acceptable, the tool should 
be more than sufRcient to state the existence of bound 
states of these systems and their approximate energy. We 
are reassured in this task by recent works dealing with 
analogous three hadron systems, for instance the KKN 
system in Ref. PT, where one could compare the results 
volving vector mesons and eventually the decuplet of baryonsof the EGA with those of more elaborate Faddeev calcula- 
For instance in Ref. [13] an explanation was given for the tions [T^] or variational calculations [22] • Also a study of 
Z\(5/2~)(1930) as a pA bound state. States of large spin the KNN system in Ref. [33] leads to similar results as 
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those obtained with variational methods using the same 
elementary amplitude for the KN interaction |^. 

The fact that we start with the pp system in 5' = 2 al- 
lows us to get bound states of the two systems studied with 
spins S = 3/2, 5/2 for the ppN and S = 1/2, 3/2, 5/2, 7/2 
for ppA. Such high spin states are difhcult to obtain with 
the right energy in quark models, which make the struc- 
ture suggested in this work a more plausible ansatz for 
the nature of these states. This is the case of the two 
Z\(5/2+) around 1740 MeV and 2200 MeV, which have 
been discussed in Ref. [TH. 



2 the interaction of vector meson and octet 
baryon 

In the Fixed Center Approximation (FC A) , we shall need 
the pN interaction and the form factor of the /2(1270) 
resonance. First we show the result of the pN interac- 
tion we used, which has been obtained in Ref. [8]. We use 
the Bethe-Salpeter equation in coupled channels (piV, wiV, 
(t)N , K*A and K*S), and the scattering matrix is given 
in matrix form as 



t = V + VGt= [l^VG]-^V, 



(1) 



where the kernel V and the meson-baryon two-body prop- 
agator G are also given in Ref. [81. In dimensional regu- 
larization the propagator G is given by 
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which is identical to that in Ref. [S], although written in 
a different form. The magnitudes mi and Mi denote the 
masses of the two propagating particles. In the present 
case they are vector mesons and octet baryons, respec- 
tively. The parameters p, and a^ are fitted with the exper- 
imental data of vector meson-baryon scattering. Here we 
shall use the values p = 630 MeV and ai(/i) — —2.0, which 
have been found in Ref. ^ to lead to a good description of 
a few vector baryon resonances (note that p and a.i(//) are 
not independent parameters and there is only one degree 
of freedom) . 



where /^r — 93A/el/ is the pion decay constant, e, e' the 
polarization vectors of the initial and final vector and /c*^, 
k '^ the energies of the initial and final vector mesons. The 
coefficients Cy are given in Ref. [8] for the different transi- 
tions between the coupled channels. This potential stems 
from the exchange of a vector meson between the external 
vectors and the baryon, involving a three vector vertex 
which is provided by the local hidden gauge Lagrangian 

m- 

For the pp system we do not need the interaction but 
only the wave function. It is obtained using the interaction 
of Ref. [25] with a four- vector contact term and diagrams 
with vector exchange in the t and u channel involving two 
three- vector vertices [12], constructing the wave function 
following the approach of Ref. [26] . 

Since the vector mesons, particularly the p and the 
K* , are rather broad, one has to take into account their 
widths. To do this, we follow Ref. [5] and replace the G 
function appearing in Eq. ^ by 



G{s, m, M ) = ^^ / dm^ x 
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where the normalization factor Nm is 
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(5) 



and the "width function" r{fri) is defined to be 



2 3 

r(m) = ro^^0(TO-mi-m2), (6) 

m^ g^n 

where Fq is the decay width of the vector meson, and 
mi, m2 are the masses of the two pseudoscalar mesons, to 
which the vector mesons decay. For the p meson, they are 
mi = TO2 = m^, and 



<?off = 



X^/'^{m'^,m1,ml) 



2m 



<?on = 



X^/^Ml^ml^ml) 
2M„ 



(7) 



For K* mesons, they are mi = m-^ and m2 = m,K, and 
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(8) 



where A is the Kallen function. 

In principle, since one has coupled channels to the piV, 

wiV, (/>A^, K*A and K*S, the three body system should 

-.-,,, . . p 1 • ,1 ,1 contain pa; iV, p0iV etc. states, in particular the pwiV state 

Under the approximation of neglectmg the three-momentum^j^g^ .^ ^^^^^^ ^^ ^^^ ^^^ ^^^ ^^^^^ pj^^ ^ ^^^^_ ^^^ 

strategy of considering all coupling channels to evaluate 
the most important amplitude and then neglect the less 
important channels in the three body system has been fol- 
lowed in similar systems, like the KNN, where the KN 



of the vector meson versus its mass, as done in Ref. [8] and 
followed here, the potential V is given by 



V^, 



1 



'4^C«.(fc 



fc")e 



(3) 
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amplitude is evaluated within the KN and irS coupling 
channels, but the nSN channels are ignored in the three 
body system (beyond what is implicitly accounted for in 
the KN amplitude) 24 . In Ref. [23 it was, however, 
shown that the explicit consideration of the irSN chan- 
nel in the three body system led to minor changes, once 
the effect of the ttS channel was incorporated in the KN 
amplitude. In the present case there are strong reasons 
to neglect the pujN, ujloN channels as explicit channels in 
the three body system. Indeed, the ujN amplitude at ener- 
gies close to the pN threshold is negligible [8, and the low 
channel couples very weakly to pp through the /2(1270) 
resonance [27] (the puj channel appears in / = 1 and is 
only relevant at energies around 1600 MeV [17]). 



3 A^-/2(1270) interaction 

The essence of the FCA is that from the three body system 
one has reasons to consider that two particles are strongly 
bound making a cluster and the third particle interacts 
moderately with the components of the cluster such as 
not to modify its wave function. Then one ignores the 
coupled channel dynamics and considers the interaction of 
the third particle with the cluster, including all multiple 
scattering steps with the components of the cluster. 

The formalism to evaluate the interaction for the N pp 
system is similar to the one used in Refs. [^l^fTl] . and 
we shall follow the same procedures and study the case 
of TV — (p/3)/2(i27o); in particular [Hll^ where one deals 
with different particles. For the N pp system, we will con- 
sider that two of the p mesons are clusterized forming 
an /2(1270) resonance, given the strong binding of the 
/2(1270) system. This allows us to use the fixed center 
approximation (FCA) to the Faddeev equations. 

The FCA to Faddeev equations is depicted diagram- 
matically in Fig. [TJ The external particle, the nucleon 
in this case, interacts successively with the other two p 
mesons which form the pp cluster. The FCA equations 
are written in terms of two partition functions Ti and T2 , 
which sum up to the total scattering matrix 

Ti = ti + tiGoT2 , 

T2 = t2 + t2GoTi , 

Tn f2{1270) =Ti+T2, (9) 



where T, 



Af/2(1270) 



is the total scattering amplitude we are 
looking for; Ti accounts for all the diagrams starting with 
the interaction of the external particle with particle i of 
the compound system; ti represents the Np unitarized 
scattering amplitude of a proton with any of the p in the 
I = pp system; Go is the loop function for the particle 
N propagating inside the compound system which will be 
discussed later on (see Eq. ([Ml). The schematic represen- 
tation is depicted in Fig. [TJ Fig. [T^) represents the single- 
scattering contribution and Fig.[TJ)) the double-scattering. 
The contributions of Fig. [T^) and b) are the two first con- 
tributions of the Faddeev equations. 





Fig. 1. Diagrammatic representation of the fixed center ap- 
proximation to the Faddeev equations for the interaction of a 
particle A with a particle B made of a cluster of two particles 
fei and 62 • Diagrams a) and b) represent the single and double 
scattering contributions respectively. In addition to these dia- 
grams there are the equivalent ones where the external particle 
scatters first with the second particle in the clusters. 

Since there are two p mesons in the present case, we 
have Ti —T2, and thus we have the following equations: 



Ti — ti + tiGoTi . 



T, 



JV/2(1270) 



3.1 Wave function of /2(1270) 



2Ti. 



(10) 



We assume that the two p mesons are bound forming the 
/2(1270), whose wave function and the corresponding form 
factor have been discussed in Refs. [^129) . It has been 
shown that the use of a separable potential in momentum 
space of the type 



V^ve{A~q)6{A-q'), 



(11) 



where A is the cutoff used in the theory for the scattering 
of two p mesons as in Ref. 116!, and q, q' are the modulus 
of the momenta, leads to the same on shell prescription for 
the scattering matrix as is used in the chiral unitary ap- 
proach. This allows to calculate form factors consistently 
with the findings of the chiral unitary approach for the 
resonances which are dynamically generated [501126] . Ac- 
cording to these references, we use the following form fac- 
tor 



FfM = 



1 



<A 
p-q\<A 

1 



(fp X 



(12) 



Mf, - 2wp(p) Mf, 
where the normalization factor J\f is 

f £p 1 

lp<A {Mf, - 2up{p)Y 



2u} p{p - q) ' 



AA 



(13) 
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of the ij system and k ^ i.j (the spectator p). This leads, 
in terms of the / — 3/2 and / — 1/2 unitarized ampli- 
tudes {rj^~ ' ' , r^~ ' '), to the following amplitude for 
the single scattering contribution: 

h-\{^t^t;"\s')^tr;"\s')), (17) 

where s' is the invariant mass of the external particle 
N and the first/second particle inside the bound state 
/2(1270): 



Fig. 2. Form factor of the /2(1270) resonance with A — 
875MeV. 



Here Mf^ is the mass of the bound state /2(1270); u!p{p) 
is the energy of the internal p meson of /2( 1270). In Fig. [2] 
we show this form factor F{q) as a function of q, where 
the cutoff parameter A is chosen to be 875 MeV. 



3.2 Single and Double Scattering Contribution 

We write the amplitude ti entering Eq. (ITUl) in terms of 
the / ~ 3/2 and / = 1/2 unitarized amplitudes (t)^ , 

iJY-~ ) . For this we consider a cluster of two p mesons 
in isospin / = 0, the constituents of which we call mesons 
1 and 2. The nucleon will be labeled with number 3. The 
two p mesons forming the /2 are in an / = state 



s' = si = S2 ^7;{s + 2TOp + Mjf - Ml 



(18) 



Now the 5'-matrix for the single scattering term is 
given by 



s\ 



(1) 



-iti 



1 1 



1 



M 



N 



M, 



N 



- V2 ^/2[j;r, V^\l EN{k) Y EN{k') 
x(27r)U{k + Kf,-k'-K'fJ, (19) 

where V is the volume of the box where we normalize 
the states to unity. The superindex of S indicates single 
scattering and the subindex that the interaction occurs 
on particle 1 of the cluster. The final expression for the 
5'-matrix for the double scattering process is 
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N 



1 



1 



1 



1 



\PP) 



1 



7=0 



Vs 



\p+p- + p-p+ + p'p") 



(14) 



EN{k) Y ENik') y/2ujp^ y2w^ y^^uj^ \/2w^ 

1 



, d^q Mn 



-hu 



U'- 



1(1,-1)) + |(-1,1))- 1(0,0))) 



where the kets in the last member of the equation indicate 
the Iz components of the 1 and 2 particles, |(/1 , /I ')). 
We take the nucleon in the state |(/1 ')) 

\P) 



(1) 



(15) 



(2^)3 '''^'EN{q) q^-EN{q)+ie 
xi27:rS{k + Kf,-k'-K'jjl^. (20) 

The subindex indicates that the first interaction is in par- 
ticle 1 of the cluster. In Eq. (PH)) . we have also taken into 
account that (fe + fc')/2 = on average. 



3.3 General form of the 5-matrix 

The general form of 5'-matrix for the nucleon-/2(1270) 
interaction can be written as 



The scattering amplitude ii in terms of the two-body 
isospin amplitudes is given by 



Si = Sft - iTNf2 (s) 
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where i and / indicate the initial and final states, respec- 
tively. By comparing the summation of S^^' and S'"^ with 
the scattering matrix S in Eq. (gl]), the iV- /2(1270) am- 
plitude can be obtained as 



TNf,{s)^2- 



yP^V^ 



1(24 



/2ujp ^2ujp> 



ti + 
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'2ujp ^J2L 
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where the notation followed in the last term for the states 
is ((/total jtotai^ I^)\tij I), where /total nieans the total isospin 



<Pq p , .Mn 



{2^f^'^'^^'EN{q^) q"-EN{q')+te'J 
2{t[+t[Go{s)t[), (22) 
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where 

Go(s) = 



2w 
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2uJn 
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with EN{q) = y^q ^ + M'^. In the center of mass system 
of the nucleon and the bound state of /2(1270), the energy 
q'^ carried by the nucleon is 



q' = 



Mj, - Ml 



5 Extension to the case of the App system 

We can follow the same procedures to study the case of 
the App system. There are just a few points we need to 
pay attention to: 

(23) 1. We need now the pA interaction. This one is obtained 
in a similar way to the pN and is done in Ref. [9]. 
There the coupled channels Ap and S*K* are taken 
into account, but the results are basically the same as 
those obtained using the single channel Ap in Ref. [T3] . 
Once again the interaction of Ap — >■ Ap stems from 

|-ig' p exchange between the two external p and the A, 
involving the three vector vertex of the local hidden 
gauge approach [25) . 
2. We need to take into account the decay widths of de- 
cuplet baryons. This is similar to what we have done 
in Eq. @. To do this, we follow Ref. |8j, and replace 

the G-function appearing in Eq. dU by G: 



2^ 



If the whole series of the expansion of Tjv/a in the fixed 
center approximation is taken into account, the scattering 
amplitude Tpff^ can be deduced from Eq. ([TU]) and we have 



T, 



JV/2(s) = Y 



2t[ 



Got[ t'-i(s')i-Go(s) 



(24) 



4 Consideration of the /2(1270) width 

The consideration of the pp — > tttt decay in L = 2 in 
Ref. [T^ gave rise to a width of the /2(1270), which has not 
been taken into account in our previous analysis. To add 
its contribution, we modify the scattering matrix T^f^ (s) 
to 
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where the normalization factor Nf^ is 
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where the normalization factor Nm is 

|.Af+2rj, 

' M-2rM 



Nm = 
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dM I I Im- 



M 



■ r'{M) 



(29) 



and the "width function" r'{M) is defined to be 
'Ai/2(M2^Mi2,M|)2Af\3 



f'(M).ro(- 



ai/2(M2,M2,M|)2A'// ^ 12;, 

(30) 
where Fq is the decay width of the decuplet baryon; Mi 
and M2 are the masses of the octet baryon and pseu- 
doscalar meson to which the decuplet baryon decays. 
Take the A baryon as an example, they are Afi = Afjv 
and Af2 = m^T^. 
3. The amplitude for the single-scattering contribution is: 

h - \ (3/2ti=;/^ + 1%"' + i/2ti=;/^) . (31) 

Besides this, we also need to change the relevant masses, 
etc. But all the procedures are similar and straightforward, 
and so we shall not discuss it further. 



and Fi{rh) is proportional to the decay width of the /2 



F.^ 



/2(1270) 



= 157 MeV, 



A(™) = r^ 



/2(1270) 



A-^/^(777^, 777^, 7773)2777/2 



V Ai 



Ai/2(,7|2 ^ rnl,m2)2m 
where the exponent in the momenta is now 2i + 1 = 5 



(27) 



6 Spin of the states 

As to the spin of systems, it is easy to determine. The clus- 
ter of pp is the /2(1270) which has J = S* = 2. Then the 
nucleon interacts with the p meson. At low energies, where 
the three-momentum of the p meson can be neglected ver- 
sus its mass, the pN — >■ pN amplitude has the form Aee' , 
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Fig. 3. Modulus squared of the unitarized N — p — p ampli- 
tude with A = 875 MeV. The solid line denotes the case with 
the /2(1270) decay width, and the dashed line is for the case 
without the /2(1270) decay width. 



Fig. 4. Modulus squared of the unitarized A — p — p ampli- 
tude with A = 875 MeV. The solid line denotes the case with 
the /2(1270) decay width, and the dashed line is for the case 
without the /2(1270) decay width. 



with e, e' the initial and final p polarization vectors [8]. 
With no extra nucleon spin dependence and L = 0, this in- 
teraction leads to degenerate pN states with spin 1/2 and 
3/2. For the case of the pA — >■ pA interaction, the same ar- 
gument holds, but now the degenerate pA states will have 
spins 1/2, 3/2 and 5/2. If we start from pp in S — 2, the fi- 
nal spin states of the three body system, which will appear 
degenerate in the approach will have J^ — 3/2+, 5/2+ for 
the ppN and J^ = l/2+,3/2+, 5/2+, 7/2+ for the ppA. 
The largest spin available would in principle be the most 
difficult to accommodate in quark models, and be better 
candidates for molecular states. 




Fig. 5. Diagrammatic representation of the formation of a 
resonance R made of a bound App and decaying into particles 
c and d. 



7 Results 



7.2 App system 



7.1 Npp system 

In Fig. [3J we show the results of |Tp for the Npp system 
with A = 875 MeV, which is suited to obtain the /2(1270) 
resonance according to the work of Ref . [TB] . We observe 
a peak around 2227 MeV with a width of 100 MeV. The 
peak does not have a standard Breit-Wigner form. The 
sharp peak could be indicative of a cusp effect but the 
threshold for ^/2(1270) is at 2210 MeV, about 17 MeV 
below the peak in Fig. [3] In fact, a small cusp peak at 
threshold is also visible in the figure at this energy. On the 
other hand, when the convolution for the mass distribution 
of the /2(1270) due to its width is considered, the peak 
of the cusp disappears, but a peak in |Tp, sfightly shifted 
to higher energies, still remains with a similar or slightly 
larger width. 

The nonstandard shape of the peak, but particularly 
the small strength of \T\^ compared to the ppA case (see 
Fig. H]) , should warn us about identification of this peak 
with a resonance. With such a small strength, other com- 
ponents of the wave function as pK* A, etc., could play 
some role. We can only take this signal as a weak indica- 
tion of a possible state around this energy in which the 
ppN component might play a relevant role. 



In Fig. 01 we show the results of jTp for the App system 
with A = 875 MeV. We observe a peak around 2372 MeV. 
Two features can be observed in Fig. |4l The peak has now 
an approximate Breit-Wigner shape and the strength of 
|rp at the peak is about 200 times bigger than in Fig. |3] 
for the N pp state, although the peak is now narrower than 
in Fig. [21 Yet, the integrated strength of the peak is still 
about 40 times bigger. The large value of T in the App 
case indicates that in a production of the resonance in 
one reaction, the magnitude of the resonance excitation 
would be large through the consideration of the interme- 
diate App state and its coupling to the resonance (see 
Fig. [5]). The consideration of the width of the /2(1270) re- 
duces the strength of the peak and increases the width of 
the resonance. The width is still relatively small, about 25 
MeV. In our approach, once the convolution for the width 
of the /2(1270) is done, the main decay channel would 
be At:!:. This is interesting to know from the experimen- 
tal point of view. This kind of information is useful since 
there is an increasing interest in looking at multimeson 
final states from the experimental point of view, and de- 
voted programs are carried out at Hall D of Jefferson Lab 
[S], at COMPASS 121] and at LEPS in Spring8/0saka, 
where Ap final states are being investigated |33) . 
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8 Discussion and conclusions 

The region of energies around 2300-2400 MeV is poorly 
known concerning N and A resonances. Although several 
resonances are compiled in the PDG [34ij it seems clear 
that there is room for many more. Concerning the N* res- 
onance obtained with 3/2+, 5/2+, there is no candidate 
in the states tabulated in the PDG. Our support for this 
state was also weak. On the other hand, the A state ob- 
tained with 1/2+, 3/2+, 5/2+, 7/2+ around 2372 MeV 
could find a counterpart in the Zi(2390)(7/2+) [31], which 
has been observed in two independent experiments [3S] 
with a mass of 2350 ± 100 MeV and in Ref. [36] with a 
mass of 2425 ± 60 MeV. It is also worth mentioning that 
this state is not seen in the latest GWU analysis [37]. 

In any case the state claimed in both 35,36 has a 
width of around 300 ± 100 MeV, which seems not com- 
patible with the narrow width of about 25 MeV of the 
state that we obtain. It is also worth recalling that masses 
of N* and A* resonances in this energy range and with 
these quantum numbers can be reached with quark mod- 
els [38) . hence maybe the narrow width is indicative of 
another kind of state, mostly made from a bound App as 
we have found. The two experiments done to claim the 
Z\(2390)(7/2+) are done with analysis of the ttN -^ wN 
reaction. As we have mentioned, the ttN is not the ex- 
pected decay channel of the resonance found, but mostly 
Atttt. a different kind of analysis would have to be done 
to eventually find this resonance. The advent of programs 
looking for multimeson final states [51113^133] offers unique 
opportunities to look for these states. The theoretical cal- 
culations, indicating the most favorable decay channels, 
can serve as a guideline for these experiments. 
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